acta physica slovaca vol. 47 No. 3/4, 1 10 



July/August 1997 



THE QUANTUM TOMOGRAPHIC ROULETTE WHEEL 1 



G. Mauro D'Ariano 2 and Matteo G. A. Paris 3 

Dipartimento di Fisica 'Alessandro Volta' dell' Universita degli studi di Pa via 
Istituto Nazionale di Fisica della Materia - Sezione di Pavia 
via A. Bassi 6, 1-27100 Pavia, Italy 

Random-phase homodyne tomography of the field intensity is a concrete example 
of the Quantum Roulette of Helstrom. In this paper we give the explicit POM 
of such measurement and compare it with direct photodetection and heterodyne 
detection. Effects of nonunit quantum efficiency are also considered. Naimark 
extensions for the roulette POM are analyzed and its experimental realization is 
discussed. 



I. Introduction 

In recent years much attention has been devoted to state reconstruction techniques, 
namely to measurement schemes providing the elements of the density matrix in some 
representation [1-4] . The problem has a fundamental interest and also practical appli- 
cations as, for example, in determining coherence properties of the field and in charac- 
terizing effective quantum interaction Hamiltonians of optical media. Nevertheless, in 
many practical situations the relevant information one aims to gain about the quantum 
state of a light beam regards the field intensity, namely the photon number. Actually, 
when looking at realistic situations, methods for precise measurements of photon inten- 
sity are especially welcome. Measurements on low-excited highly-nonclassical quantum 
states of radiation represent, in fact, a difficult task, due to limitations of currently avail- 
able photodetectors. On one hand, it is very difficult to discriminate single photons in 
any desidered range of intensity. On the other hand, in the low photon-number region 
(n < 10) one can discriminate a single photon, however with low quantum efficiency. 
Therefore, alternative methods to detect intensity need to be considered, whereas their 
performances and their robustness to low quantum efficiency should be compared in the 
very-quantum regime, versus the limitations of customary photodetection. 

A way to circumvent photodetection problems is that of amplifying the the output 
photocurrent by mixing the state under examination with a highly excited coherent 
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state. This is the basis of homodyne and heterodyne detection schemes, and, in turn 
homodyne detection is the basis of tomographic state-reconstruction. In this paper we 
address the measurements of the field intensity by random-phase homodyne tomography 
and compare it with heterodyne detection. 

Random-phase homodyne tomography has been suggested [6] as an effective tool to 
measure photon intensity, and after has received some theoretical attention [5]. More 
recently experiments have nicely demonstrated its ability in determining photon statis- 
tics for highly nonclassical states of radiation [7]. Here, we derive the probability op- 
erator measure (POM) that describes the random-phase homodyning of the intensity, 
regarding this measurement as a concrete example of the quantum roulette of Helstrom 
[8]. This abstract framework is suitable to discuss different Naimark extensions, thus 
suggesting novel experimental realizations. 

The paper is structured as follows. In Section II we describe the measurement 
of the field intensity by the quantum tomographic roulette. We derive its POM and 
discuss its precision. In Section III the tomographic roulette is compared with direct 
photodetection and heterodyne detection also in presence of nonunit quantum efficiency. 
In Section IV we analyze the possible Naimark extensions to the roulette POM, whereas 
Section V closes the paper with some concluding remarks. 

II. Measuring Intensity by Quantum Roulette 

The concept of Quantum Roulette has been introduced by Helstrom in his book [8] 
about twenty years ago. A Quantum Roulette is described by a POM of the form 

M 

ft m = £z fc £W m = l,...,N, (1) 
k=i 

where Zk > 0, J2k=i Zk = -0 an< ^ tnc ^«?' s are families of orthogonal projectors corre- 
sponding to different observables labelled by k, say 0^ k \ in formula 

N 

Eg>EW=5 mn EW £4? = L (2) 

m— 1 

Each experimental event corresponds to the measurement of one of the observables 
0^ k \ chosen at random according to the probability distribution Zk- At the time of 
Hclstrom's proposal the Quantum Roulette was just a tool to illustrate an example of 
generalized measurements that do not correspond to selfadjoint operators. Nowadays, 
such measurement can be now performed in quantum optics labs. Let us consider the 
homodyne detection of a nearly single-mode radiation field. When the phase <j> of the 
local oscillator is fixed, the field-quadrature x$ = ^(aJe 1 ^ + ae~ l< ^) is detected and the 
measurement is described by the POM 

dE^\x) = \x) u {x\dx, (3) 

\x)<p denoting eigenstates of x^. The tomographic detection of the field corresponds 
to scan the local oscillator phase over [0, n]. When the phase <j> of the local oscillator 
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is unknown one deals with a Quantum Roulette, each experimental event being the 
measurement dE^ (x) with random <f>. The experimental outcomes are distributed 
over the whole real axis according to the probability distribution p{x) = Tr{gdjl(x)} 
where dfi(x) is the nonorthogonal POM 

dfi(x)=dxf — |x)00(x| . (4) 
Jo * 

Inserting the number state expansion of \x)<f, 
in Eq. (4) we obtain the POM of the roulette 

™-^±^™-*&-~^ 

n— v ' 

where H n (x) denote Hcrmite polynomials. 

The POM dfl(x) is an operator function of the number operator only and the outcome x 
will be an estimate-generally biased-of the field intensity. We now proceed in deriving 
an unbiased estimate. As it was shown by Richter [9] the expectation value of any 
normally ordered product (a tn a m ) can be obtained from tomographic data by averaging 
the kernel integral 

1l[<J n a m ]{x,<j)) = e <*(">-») H n+ (V2x) 



2(n+m)/2 ' 71 + 771 



(7) 

m J 

over the probability distribution p(x, <j>). For the number operator h Eq. (7) defines the 
kernel 

y = TZ[a^a}(x) = 2x 2 -^, (8) 

which is a phase-independent quantity, hence is suitable for estimation by a Quan- 
tum Roulette. The quantity y traces the field intensity by averaging over the roulette 
outcomes distribution 

/oo 
dxTZ[a)a]{x)p(x) = (h) . (9) 
-oo 

Indeed, Eq. (9) shows that th function y(x) in Eq. (8) is the unbiased field intensity 
estimator for Quantum Roulette. The single outcomes y = lZ{a^a](x) are random 
numbers distributed over the interval [— 1/2, oo]. It is clear that the determination in 
Eq. (9) is meaningful only when also a its statistical deviation is specified. The latter 
is given by 



Ay = vV-y 2 ' ( 10 ) 
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where 



y 



dx 7£ 2 [ala](x) p(x) . 



The explicit expression of the statistical deviation is given by 

(n 2 ) + (h) + ij , 



Ay 2 = (An 2 ) + 1 



(11) 



(12) 



(An 2 ) being the intrinsic photon number fluctuations. We can also specify the whole 
probability distribution p(y). In fact, from the Radon-Nikodym derivative of the 
roulette POM in Eq. (6) we arrive at the roulette POM for field intensity 



dfi{y) = 



dy e-fa+V 2 ) ifl,(Vy+l72) 
V* y/yTT/2 2 ata (ata)! 



(13) 



III. Quantum Roulette versus Photodetection and Heterodyning 



As it emerges from Eq. (12) the Quantum Roulette measurement of the field intensity is 
noisy, as compared with ideal photodetection. Here we analyze its performances in the 
realistic case of nonunit quantum efficiency. When dealing with r\ < 1 the overall output 
noise (An 2 ),, is larger than the intrinsic quantum fluctuations (An 2 ) which represents 
the minimum attainable noise in a measurement of the intensity. For non unit quantum 
efficiency, the roulette POM becomes a Gaussian convolution of the POM (6) 



dy' 



27TO- 2 



dfi(y') exp 



'\2 



(y-y' ) 

2a 2 



1 — 7] 

4:7] 



whereas the unbiased estimator is now given by [10] 

y r , = n v [a)a]{x) = 2x 2 - ^~ 

Zrj 



Inserting Eqs. (14) and (15) in Eq. (10) one has 
A^ 2 = (A^) + i(n^) + (n) (- 



2 3 
V 2 



1 

v 



(14) 



(15) 



(16) 



This noise has to be compared with the rms variance of direct detection for nonunit 
quantum efficiency, which is given by 



(An 2 ), = (An 2 ) + (n) 



1 



1 



(17) 



The difference between Ay n and {An 2 ) v defines the noise N R [n] added by Quantum 
Roulette with respect to direct detection 



N R [n] 



1 



(n 2 ) + (n) 



1 + 



1 



(18) 
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Nn[n] is always positive, and roulette determination is more noisy than direct photodc- 
tection even for nonunit quantum efficiency. 

Let us now consider heterodyne detection, namely the joint measurement of two com- 
muting photocurrents, which, in turn, trace a pair of conjugated field quadratures. 
Each experimental event corresponds to a point in the complex plane of the field ampli- 
tude and these outcomes are distributed according to the generalized Wigner function 
W s (a, a) with ordering parameter s related to the quantum efficiency as s = 1 — 2ry~ 1 . 
Starting from the relation 

lap= / d 2 aaa* WJa,a) = (a>a) + - , (19) 
Jc V 

we are led to consider the shifted square modulus I v = \a\ 2 — 1/n as the heterodyne 
unbiased estimator for the field intensity. For unit quantum efficiency r\ — 1 heterodyne 
detection measures the Husimi Q-function (a|£|a), and thus the POM for the field 
intensity dfi(I) is the marginal one of the Arthurs-Kelly coherent-state POM dfi(a) — 
7r _1 |a)(a|, we have 

dfi(I) =dlY e-( /+1 ) iL±2l. \ k ) = dl e -< /+1 ) (/ + 1} ° ta . (20) 
^ kl (ato)! 

For the determination I v = (h) we need to specify the statistical deviation 

^ = P r -l q \ (21) 
By using Eq. (19) and the following relation 

H* = f d 2 a a 2 a* 2 W s (a, a) = {rt) +(--l)(n) + ±, (22) 

we arrive at the result 

A7f - (A^) + (23 ) 
\V J V 

which represents the precision of heterodyne detection in measuring the field intensity. 
From Eqs. (17) and (23) we obtain the noise added by heterodyne detection with 
respect to direct detection 

(24) 

-/Vy[n] is always a positive quantity, thus also heterodyne detection is more noisy than 
direct detection for any value of the quantum efficiency. 

A direct comparison between Quantum Roulette and heterodyne detection can be ob- 

2 2 

tained by considering the difference A RH [n] = Ay n — AI V . From Eqs. (16) and (23) 



N H [n] = - 



(n) + - 
V 
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one has 



A RH [n] = (A^)_(A/2} = 



(«=») - (n) - 



n=0 



(25) 



being the diagonal elements of the signal density matrix. A^#[«] has no definite 
sign, when changing the states of radiation. Therefore, it is matter of convenience 
to choose between the two kinds of detection scheme, depending on the state under 
examination. For any value of the quantum efficiency rj the quantity [n 2 — n — 
becomes positive for n larger than the threshold value 



riT 



1 + y/1 + 4/r/ 2 



(26) 
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Comparison between the Quantum Roulette and heterodyne detection noises in the determi- 
nation of the field- intensity. The lines where A RH [h] = for squeezed states are plotted for 
different values of the quantum efficiency (from left to right: rj — 1.0,0.75,0.5,0.25 and 0.1) 
as a function the total mean photon number iV and the squeezing photon fraction /?. The 
region on the left of each curve corresponds to states for which the Quantum Roulette is more 
convenient than heterodyne detection. This region becomes larger for decreasing rj. 

This means that the Quantum Roulette is more precise than heterodyne detection 
for low excited states, as the higher-n terms get a lower weight g nn . Moreover, the lower 
is the quantum efficiency the larger is the threshold value, and hence the region where 
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the Quantum Roulette is convenient. This means that the Quantum Roulette is very 
robust to low quantum efficiency. For coherent states the Quantum Roulette becomes 
convenient with respect to heterodyne detection when the field intensity is lower than 
the value (n) = r\~ x . For squeezed states (we consider the case of zero signal and zero 
squeezing phases) Eq. (25) can be rewritten as 

A RH [n] =N 2 + 2pN{l + 0N) + (1 - 0)N(1 + 2[3N + 2y/0N{l+ 0N)) -N-\, 

where TV denotes the total mean photon number of the state and (3 the squeezing photon 
fraction, namely the ratio between the (mean) photons engaged in squeezing and the 
total one (for j3 = we have a coherent state and for (3 = 1 we have squeezed vacuum). 
In Fig. the lines where Arm [n] is zero are plotted for different values of the quantum 
efficiency. The region on the left of each curve corresponds to states for which the 
Quantum Roulette is more convenient than heterodyne detection. 

IV. Naimark Extensions 

Naimark theorem [11] assures that every POM is a partial trace of a customary or- 
thogonal projection- valued measure on a larger Hilbcrt H = Hs ® Tip space which, 
itself, represents the original system interacting with appropriate probe mode(s). This 
extension is not unique, corresponding to the different possible physical implementation 
of a given measurement. 

For the abstract Quantum Roulette of Helstrom, defined in Eq. (1), there exists a stan- 
dard recipe to obtain a Naimark extension. This is found by considering the following 
projectors [8] 

M 

Em^Y. ^ ® hbXwfcl, (27) 
k=i 

where \u>k) denotes an orthogonal set of states in the extension (probe) Hp Hilbcrt 
space. In fact, by preparing the probe in the state 

M 

hM = £4 /2 K> ; (28) 
k=i 

we have 

Trpjis® \^p){ME m ) =il m - (29) 

The POM in Eq. (4) is the continuous version of the quantum roulette with <p playing 
the role of the label k. The Helstrom recipe to achieve a Naimark extension requires 
an orthogonal POM for the phase (f> in Hp yielding a resolution of identity in [—n, n]. 
As it describes a phase variable, such a POM has to satisfy the additional requirement 
of covariance, namely 



27T 



(30) 
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S being the phase-shift generator in the probe Hilbert space, and P a suitable positive 
operator. 

The "minimal" implementation for the roulette POM would be a two-mode system, 
with the additional mode playing the role of the probe. However, such an imple- 
mentation cannot be fully quantum, because no single-mode orthogonal set the phase 
variable exists. By using the canonical (Susskind-Glogower) phase POM dfi((f>) = 
{Zn) d<f>^2 nm exp [i(n — m)4>] \n)(m\ we have an approximated extension described by 
the two-mode (non-orthogonal) POM 

dM(x) — dx [ \x) l p r p(x\ <g> exp [i(n — m)(f>] \n)(m\ , (31) 

J 27T ' 

nm 

which corresponds to the measured photocurrent 

X = j x dM{x) = a f e_ + ae+ , (32) 

e± being the raising and lowering operators on Hp 

e+ = fet 1 g_ = - p b . (33) 

An exact two-mode implementation can be achieved in the semiclassical limit of highly 
excited probe mode. The basic idea is that coherent states \z) provide exact phase 
states for high amplitude \z\ — * oo. In this limit we can neglect the ordering in (33) 
and writing 

a 1 e_ + ae + ~ : ~ , — , (34) 

+ Vtfb + i y/jpwj v ; 

which coincides with the homodyne photocurrent. The prescription (28) for the probe 
preparation becomes 

g P = hm f f- ||z|e^)(|z|e^| . (35) 

Eqs. (34) and (35) unambiguously identify random-phase homodyne detection as the 
minimal (semiclassical) implementation of the Quantum Roulette POM of Eq. (4). 
A fully quantum extension requires more than one mode as the probe. An example is 
provided by the heterodyne phase eigenstate, which are defined on a two-mode Hilbert 
space [12]. 



V. Conclusions 

The Quantum Roulette is a concrete example of the Quantum Roulette of Helstrom. In 
this paper we compared intensity measurement by the Quantum Roulette with that from 
direct photodetection and from heterodyne detection. Direct photodetection, though 
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only in principle, remains the most precise detection scheme, also in the case of nonunit 
quantum efficiency. On the other hand, the choice between Quantum Roulette and 
heterodyne detection is a matter of convenience, depending on the state under exami- 
nation. For coherent states, the Quantum Roulette becomes convenient with respect to 
heterodyne detection when the field intensity is lower than threshold value (n) = n^ 1 . 
In general, Quantum Roulette is more convenient for low excited states, and is more 
robust in presence of low efficient photodetectors. 
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